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Abstract 



Heavy baryon chiral perturbation theory is extended to include the ef- 
fects of quenching. In this framework the leading nonanalytic dependence 
of the heavy baryon masses on the light quark masses is studied. The size 
of quenching effects is estimated by comparing the results of quenched and 
ordinary chiral perturbation theories. It is found that in general they can 
be large. This estimate is relevant to lattice simulations of the heavy baryon 
masses. 
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Typeset using REVT^ 



I. INTRODUCTION 



Recent technological developments have made it possible to perform lattice simulations 
of QCD with increased accuracy. At the same time the precision of the lattice calculations 
depends on the understanding of the errors due to various approximations used in lattice 
simulations. One of the most important errors is due to quenching. The modifications 
induced by quenching at short distances have been studied in systems of quarkonia 
and are more or less well understood. Because of its nonperturbative nature, the long 
distance effects of quenching are more difficult to quantify. One approach is to use Quenched 
Chiral Perturbation Theory (QChPT). The idea was first proposed by Sharpe 03 , with the 
formalism further developed by Bernard and Golterman . The advantage of using the 
quenched version of Chiral Perturbation Theory (ChPT) is that its predictions follow from 
the basic properties of the underling theory, and ChPT describes the low energy dynamics 
of ordinary QCD very well. 

More recently, QChPT has been extended to describe the interactions of soft pions with 
baryons heavy mesons P,p!0| and vector mesons [0]. Here it is further extended to 

describe the interactions with heavy baryons. Ordinary chiral perturbation theory for heavy 
baryons has been formulated by a number of authors |ll2|- |T^] . We adapt their formalism to 
quenched QCD. 

As every effective theory, ChPT contains a number of undetermined couplings all of 
which must be fixed from experiment. Only for the self-interactions of pions are there 
enough data to carry out this procedure beyond the leading order. The situation is even 
worse in the quenched case, where the only "data" available are extracted from lattice 
simulations. For this very reason we cannot obtain any accurate quantitative predictions 
from our calculations. Nevertheless, there are certain loop corrections whose existence is 
predicted unambiguously by the lowest order Lagrangian, in particular, terms that have 
nonanalytic dependence on the light quark masses. Contributions of these terms to the 
matrix elements cannot be absorbed into the higher order terms of the chiral expansion. 
Furthermore, since the nonanlytic dependence arises from the infrared region of the loops, 
it is particularly sensitive to the low energy behavior of the theory. 

In this paper the nonanalytic corrections to heavy baryon masses are computed in the 
framework of QChPT. Comparison with the analogous corrections obtained in ChPT shows 
the qualitative differences between the quenched and unquenched cases. It is also found 
that these nonanalytic terms could in general be large. This is important for the lattice 
calculations of the heavy baryon masses, where results are extrapolated linearly in the light 
quark masses. 

Since quenched QCD is not a unitary theory, it is possible that QChPT does not describe 
its behavior at low energies correctly. However, the most important qualitative feature of 
QChPT, the presence of quenched logarithms of the form Mplnm^, was first derived in 
strong-coupling perturbation theory [QJ^ . Moreover, the coefficients of these terms obtained 
by the two different approaches are the same. We thus expect that the description of the 
low energy limit of quenched QCD by QChPT is valid. 
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II. CHIRAL PERTURBATION THEORY FOR HEAVY BARYONS 



Heavy Quark Effective Theory (HQET) and ChPT are combined to describe the inter- 
actions of heavy baryons with soft Goldstone mesons. Here we give a brief review of the 
formahsm. 

To leading order of chiral expansion, self-interactions of soft pseudoscalar mesons are 
described by the Lagrangian 



^ (Tr[9^S9^St] + 45o Tr[M 

o ^ 

where S = and ^ = e**^/-^. Eight light pseudoscalar mesons are combined in a matrix 

76^ 



(2.1) 



V 



vr 

K- 



'V2 



71+ 


K+ \ 











(2.2) 



The normalization is such that ~ 130 MeV. The leading symmetry breaking term 
depends on 
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MH. = -(ew+eMo 
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where M is the quark mass matrix: 



M 



( fTlu 
TUd 

V 
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The mass term of ( |2.1| ) breaks the SU{3)l x SU{3)r symmetry which is restored in the 
limit muiTndiTns 0. In this limit the eight pseudoscalar mesons become massless. Note 
that the meson masses are nonlinear in the quark masses. 

Heavy baryons contain one heavy quark, c or 6, and two light quarks. According to 
HQET, in the leading order of heavy quark mass expansion, properties of the heavy baryons 
depend entirely on the configuration of light quarks and are independent of the spin and the 
flavor of the heavy quark. The three light quarks, u , rf, and s, that form the fundamental 
representation of flavor S'f/(3) can be combined in pairs to form an SUf{3) antitriplet and 
an SUf{3) sextet. There is a correlation between the spin wavefunction and the flavor 
wavefunction of the light quarks; those in the SUf{3) antitriplet combine in a spin state and 
those in the SUf{3) sextet have spin 1. The former, combined with a spin 1/2 heavy quark, 
form a triplet of baryons with = The spin 1 combination of light quarks together 
with the heavy quark forms two sextets of baryons with = |^ and = which are 
degenerate in the infinite heavy quark mass limit, Mq — > oo. Since the configuration of 
light quarks in the triplet baryons differs from that in the sextet baryons, the two are not 
degenerate even in the Mq oo limit. 

The two sextets of baryons can be described with one field S*'/ • 
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qij _ qij I qij 
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where is the 4-velocity of the baryon. The baryons with = are combined in a 



matrix Bq 
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baryons. Here the 



The matrix B^^ is similar to Bq, except that it contains J — 2 
notation of Ref. is used, where the superscripts denote the Is projection of the isospin 
of the baryons. 

The baryons in the three dimensional representation of SUf{3) are described by a field 

rpij ^ 



(2.7) 
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To lowest order in the 1/Mq expansion, strong interactions of the heavy baryons with 
the soft pions are described by the chiral Lagrangian Q 

Ain = ^ tr[T {vV)T]-t tT[S^ {v ■ V) 5^] 



Ao 



tr[TT] + Ai tr[5^M+5^] + A2 tr[5''5^] tr[M4 

+ y tr [TM+T] + y tr [TT] tr[M+] 
C,nt = 93 (tr[TA^5,J + h.c.) + tgMS'A^S''] v'e^,^, . 



Here 



is a covariant derivative, with a vector current 



(2.9) 



(2.10) 



^The normalization and the notation for coupling constants are the same as in |16|. 
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V, = l{^d,e + ^^d,^) , (2.11) 
and the axial current, A^, is defined as follows: 

A, = l{^d,e-ed,^) . (2.12) 

The covariant derivative of T*-' has the same form. 

We chose to absorb the mass of the sextet baryons into the static phase of the heavy 
baryon fields. The mass splitting between the sextet and triplet baryons, finite as Mq oo, 
is denoted by Aq. The propagators for = I"*" and = sextet baryons are 



2 vk 

and 



(2.13) 



'2 vk ' ^ ^ 



respectively, while the propagator of the triplet baryon is 

1 + «f i 



(2.15) 



III. QUENCHED QCD 



Quenched QCD is obtained from the ordinary theory by removing the disconnected 
fermion loops. Formally this is done by introducing a bosonic "ghost" partner for each 
quark g*, which has the same mass and couples to the gluons the same way that the quark 
does. The sign difference between the fermionic and the bososnic loops results in cancellation 
of the two. The symmetry of the quenched theory is enlarged from SU (3)^ x SU (3)r to the 
semi-direct product {SU{3\3)l x Sf/(3|3)R)®f/(l). 

In this larger theory the Goldstone matrix becomes a supermatrix. 



n 



X 7f 



(3.1) 



where the quark/ghost content of the fields is vr ~ qq, x'' ^ QQ, X ^ ^''^^ ^ ~ QQ- Each of 
these is an 3 X 3 matrix; for example 



TT 



TT 

K- 
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(3.2) 



Note that x and are fermionic fields, while vr and vr are bosonic. 
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Notice one important difference from ordinary ChPT: the inclusion of the SU{3) singlet 
in the theory. In ordinary QCD the anomaly pushes the mass of the rj' up to the chiral scale 
and it decouples from the low energy effective theory. In the quenched case, because of the 
absence of disconnected quark loops, this decoupling does not occur and the rj' has to be 
included in the theory. In the quenched case additional terms are required to describe the 
dynamics of the anomalous field. The lowest order Lagrangian is similar to Eq. 

^Qx = y (Str[9^S9'^St] +45o Str[^+]) 

+ ^ (Aq Str[9^n] Str[a'^n] - Str[n] Str[n]) , (3.3) 

with S = ^ = e*n/^ and 

is defined analogously to M+. The "supertrace" Str is defined with a minus sign for 
the ghost-antighost fields. Normalization of Aq and Mq is such that they have no implicit 
dependence on Nf, the number of flavors. 

1/2 ~ 

The presence of the Str[n] = A^^' {rj' — r]') term leads to a double-pole structure for the 
propagators of the flavor-neutral mesons. In the basis where these mesons correspond to 
gj^j and (jiq^ this propagator takes the form 

where ml = 2i?ow^i, and = 1 if z corresponds to a quark and ej = — 1 if z corresponds to a 
ghost. The second term in the propagator can be treated as a vertex, called a hairpin, with 
the rule that it can be inserted only once on a given meson line. 

Inclusion of quenching effects in the heavy baryon sector is straightforward. The baryon 
matrices are promoted to supermatrices and Tr is replaced by Str. The kinetic and the mass 
terms of the quenched Lagrangian of the heavy baryons are 



A^i^ = ^ Str[T iv.V)r]-^ Str[5^ (v ■ V) S, 
'C-^ls = ^ Str[Tr] + Ai Sti[S''M+S,] + A2 Str[5%] StT[M+] 

+ ^ Str[TA^+T] + ^ Str[TT] Str[;W+] . (3.6) 



where 



with Kn and describing the heavy baryons obtained by replacing one of the quarks by 
its "ghost" partner and is the heavy baryon with both light quarks replaced by "ghosts" . 
The superfield T is defined the same way. 
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T=[l f). (3.8) 
The covariant derivative, the vector and the axial currents are defined analogously to 

im) - (HI- 

In the quenched case, Str does not vanish and the symmetry allows the baryon fields 
to couple to the fiavor singlet Goldstone mesons through this term. Thus the interaction part 
of the quenched Lagrangian has one additional coupling that is absent in the unquenched 
case: 

C?,,^ = 9s (Str[T^'^5^]+h.c, 

+ igi Str[5''5'^] Str[^''] v^e^^y^, . (3.9) 



IV. NONANALYTIC CORRECTIONS TO HEAVY BARYON MASSES 

We encounter two types of integrals in our calculations. One, Jf^, arises from the chiral 
loops without the hairpin vertex and the other, from the loops with the hairpin vertex 
on Goldstone meson line: 



J (zvr) [p'^ — ni'^ + ie)[v • (p + k) + Aq + le) 



and 



jr(x, 2/, Ao,t; ■ k) = Um^ - Aox')Jr{x,Ao,v- k) - {Ml - A,y^)J^^ {y, A^.v ■ k) 



^2 _ y2 



(4.2) 



We truncate the series to include terms of order m^, mqlogm^ and m^/^ but not those 
of order and higher. If not otherwise indicated, dimensionful parameters Aq and Mq are 
considered to be of order rriq. With these assumptions, terms contributing to the baryon 
mass corrections are:0 

jr{m, 0, 0) = {vf'v'' - g^'')h{m) (4.3) 



where 



= (4.4) 

h(x, y) = [{MS ~ A„x^)h{x} - (Ml - A„;/)h(y) 

X y 



^Setting Aq = in ( [4.1| ) and ( [4.2[) means that we neglect terms of order AQiriqlii'i 
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FIG. 1. Diagrams contributing to the flavor diagonal sextet baryon mass corrections. The thick 
sohd hne denotes the sextet baryons and the dashed hne denotes the Goldstone mesons. 
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FIG. 2. Diagrams that modify the masses of sextet baryons with hght quarks of two different 
flavor. The thin solid line denotes the triplet baryons. 

I2{x,y) has the hmit 

h{x,x) = (^^M^x - ^Aox^) . (4.5) 

Here we only compute corrections to masses arising from the chiral loops and not from 
the 1/Mq expansion. Chiral corrections do not lift the degeneracy between the Bq and Bq 
baryons. However, the baryon masses inside the sextets as well as inside the triplet are split. 
Mass corrections to the baryons that have the same flavor light quarks differ from those 
that have light quarks of two different flavors. The latter get contributions from loops which 
contain triplet baryons. This is different from the unquenched case, where the loops with the 
triplet baryons contribute in both cases. The reason for this difference is that the diagrams 
with the triplet baryons which do not contain the hairpin vertex require a closed quark loop, 
and in the quenched theory they vanish. On the other hand, the hairpin vertex appears only 
on the flavor diagonal meson line and these mesons only couple the triplet baryons with the 
off-diagonal sextet baryons. The diagrams that modify sextet baryon masses are shown in 
Fig. ID and Fig. |^. The contributions from these diagrams are: 

SMfi = 2g2gih{mii) + gjhimu, niu) , (4.6) 

for the diagonal baryons, and 
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FIG. 3. Diagram that modifies the masses of triplet baryons. 

SM^j = g2gi ihimu) + hinijj)) + ^gj {himu, ma) + himjj, nijj) + 2/2(771^^, rujj)) 

+^gl ihimu, ma) + himjj, m^j) - 2/2(7/1^^, nijj)) (4.7) 

for the off-diagonal sextet baryons. hix) and hix^y) are defined in Eq. (|4.4| ). In the limit 
of exact isospin symmetry the term proportional to g1 in Eq. ( [4.7D vanishes when the baryon 
contains up and down quarks. 

As was discussed above, in the case of triplet baryons the only diagrams that contribute 
to the mass corrections are those with the hairpin vertex (Fig. |]). The triplet baryon mass 
corrections have the following form: 

3 

8Mf- = -gj ihimu, ma) + himjj, rrijj) - 2himii, rrijj)) (4.8) 

This expression has the same form as the correction to the off-diagonal sextet baryon 
masses that arise from the diagrams with the triplet baryons. Thus in the limit of exact 
isospin symmetry the mass corrections to Ag baryon vanish. 

There is no mixing between the triplet and the sextet baryons in the leading order of the 
1/Mq expansion. Mixing would require the breaking of heavy baryon spin symmetry, which 
is exact at this order of the heavy quark expansion. 

Finally, we note that these results are readily generalized to an arbitrary number of light 
flavors. 



V. PHENOMENOLOGY 

As was mentioned in the introduction, there are many undetermined parameters in the 
theory. In the quenched case there are additional parameters that are absent from the 
Lagrangian of ordinary ChPT. Of course the parameters of the quenched and unquenched 
theories, even those that are present in both cases, might not be the same. Here couplings 
that are present in ordinary ChPT are assumed to have the same value in the quenched 
case. Other couphngs we assume to be of order one. Calculations are done for two sets of 
values for Mq and Aq: Mq = 400 MeV, Aq = and Mq = 100 MeV, Aq = 0.2. For more 



detailed discussion of the values of these parameters we refer the reader to |TT|. 

There is a nontrivial relation between the sextet baryon masses that holds in both 
quenched and unquenched cases in the leading order of SUjiS) breaking. Namely, in the 
limit of exact isospin symmetry: 
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5M = + Muq - 2Ms'^ = . (5.1) 

This relation gets corrections in the next order of the chiral expansion. In the hmit of 
vanishing Aq, this correction is [|T6| 

Msg + Mn^ - 2Ms^ = ^^{gl - 9l){^rn% - SmJ - ml) , (5.2) 

assuming the tree level relation between the Goldstone mesons, Arrij^ — 3m^ — = 0. In 
the quenched case the equivalent relation takes the form 

Ms^ + - 2Ms'^ = 1(^,2 _ g^) [l2{muu, rriuu) + h{.mss, m^s) - 2/2(m„„, m^^)] • (5.3) 



To simplify the last relation let us consider the limit = = 0. Then Eq. ( ^.3|) 
becomes: 

+ - 2M=^ = -^^{gl - gl) [Mlrriss + A,ml] . (5.4) 



The main difference between the two cases is the presence of the term proportional to 
in the quenched case. The Goldstone meson masses are linear in m^^. Fig. shows 



1/2 



the dependence of these corrections on while keeping the value of the strange quark 
mass fixed. (This is the extrapolation used in lattice calculations to extract the values of 
the baryon masses.) In both cases the corrections to 5M are small but they have different 
signs. The plot indicates the very different behavior of the nonanalytic corrections in the 
two cases. 

We also investigate the mass splittings between the sextet and the triplet baryons. The 
splittings are calculated in the limit of exact isospin symmetry.^ As was discussed in the 
previous section, in this limit Aq does not receive corrections at the first order of the chiral 
expansion. The results for the mass splittings are:^ 



Aq = Ao - 2^72^1 - gi^, [-^M^rn^ - -A,mi ) (5.5) 



and 

1 



Sq - Sq = Ao + g2gi{Ii{m^) + himss)) + ^glihim^, m^) + hirUss, m^J + 2/2(m^, m^^)) 
1 



-^gl{l2{m^, m^) + hiniss, niss) - 2/2(m^, m^^)) . (5.6) 



^In this limit muu = fndd = n^ir- 

"^Not included in the following expressions are the contributions from the mass terms of the 
Lagrangian (|3.6| ) that are linear in the light quark masses. These terms depend on the undetermined 
parameters, Aj. We found that these contributions are small if one considers the linear combinations 
of Aj's to be of order 1. 
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FIG. 4. Dependence of 6M on in the quenched and unquenched cases. The sohd hne 
corresponds to the unquenched case, the dashed hne to the quenched case with Mq = 400 MeV 
and Aq = 0, and the dashed-dotted hne to the quenched case with Mq = 100 MeV and Aq = 0.2. 



The last relation simplifies in the rriu = ma 



limit: 



An- 
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92+9i) - 9291 



(5.7) 



487r/2 ^2"^ ' "^"n27r/2- 

For simplicity of the analytic expressions only the cases with = are shown. 

We would like to point out two things. First, the corrections to the mass splittings, as the 
correction to (5M, display nonanalytic dependence on the light quark masses that is different 
from ordinary ChPT. In the quenched case the leading nonanalytic term is m^^. Notice 
that the coefficient of the term linear in the Goldstone meson mass depends, apart from the 
Mq term, on the parameters that exist in the unquenched theory. One can estimate the size 
of this term assuming that the corresponding parameters in the quenched and unquenched 
theories are the same. Second, there are terms proportional to gi, the coupling that is 
absent from the ordinary theory. In general this term gives a large contribution to the 
mass corrections, which depends on the value of gi, of course. This makes the estimate of 
quenching errors even more uncertain. The dependence of the mass splittings on for the 
fixed strange quark mass is shown in Fig. |^. Here gi is taken to be 0.5 and Aq = 200 MeV. 
Note that, by contrast, the terms proportional to gi are absent from the corrections to 
5M. In this particular combination of the heavy baryon masses these terms cancel. This 
indicates that one has to be careful when using such mass relations for the estimate of 
quenching errors. 

Recently these splittings, as well as the heavy baryon masses, were calculated on the 
lattice ||l^| . In that analysis, calculations are performed for the light quark masses of order 
of the strange quark mass, and then the results are linearly extrapolated to the chiral limit of 
the up and down quark masses, while the strange quark mass is fixed to its phenomenological 
value. In particular, the following relation for the mass of the heavy baryon is assumed: 
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0.00 0.10 0.20 0.30 0.40 0.50 
m.l [GeW ) 

FIG. 5. Dependence of the sextet-triplet mass splittings on m^. The solid and dashed lines 
correspond to the E'q-Eq splitting with Mq = 400 MeV , = and with Mq = 100 MeV , Aq = 0.2 
respectively. The line with small dashing and the dashed-dotted line correspond to the Sg — Ag 
splitting with Mq = 400 MeV ,Ao = and with Mq = 100 MeV , Aq = 0.2 respectively, gi is taken 
to be 0.5. 

MQij = fx{MQ) + C{mi + mj), (5.8) 

where Mq is the mass of the heavy quark and rrii and rrij denote the masses of the hght 
quarks. ^{Mq) is the mass of the baryon in the chiral limit, = = = 0. /i and C 
are different for the sextet and triplet baryons. The mass splitting between the sextet and 
triplet baryons is: 

AMgi,- = Ao(Mq) + D{mi + m,) . (5.9) 

From Eq. (|5.9|) it follows that in the chiral limit for up and down quarks m„ = = 0, 

Ms^-Ma^=Ao, (5.10) 

and this agrees with Eq. ( p.5|) when m.^- = 0. In ChPT Aq is a free parameter and has to 
be determined from the fit to the experimental data. For our choice of Aq = 200 MeV the 
quenched result is in agreement with the lattice result. 

The situation is different for the Sg — Sg mass splitting. The quenched result is given 
by Eq. (|5.7|). The lattice result is 

- Mh^ = Ao(Mg) + Dm, (5.11) 

In this case the finite value of the strange quark mass results in the modification of the tree 
level splitting even in the limit = uid = 0. In QChPT the Sg — Sg splitting, in contrast 
to the Eg — Ag splitting, has different values for the two sets of Mq and Aq parameters in this 
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limit. As Fig. ^indicates, the corrections to the — Eq sphtting could in general be large. 
The plot also shows that the size of the corrections strongly depends on the parameters Mq 
and Aq. 

The most important result, however, is the different dependence of the lattice and 
quenched results on the light quark masses. The nonanalytic terms, present in QChPT, 
are not accounted for in lattice simulations. The presence of these nonanlytic terms is un- 
ambiguously predicted by the quenched theory and they can be calculated from the tree 
level Lagrangian. Because of the presence of several undetermined parameters, it is possible 
that the coefficients of the nonanalytic terms are small. The only way to determine the size 
of these terms is to take them into account in lattice simulations and perform the nonlinear 
extrapolations. 

VI. SUMMARY 

Let us summarize briefly the results of our calculation. The main qualitative result 
is the dependence of the heavy baryon masses on the square root of the light quark mass 
m^/^. This is different from ordinary ChPT where the leading nonanalytic dependence goes 
as m^/^. The presence of these terms is unambiguously predicted. Their coefficients are 
calculated from the tree level Lagrangian and they can not be modified by the inclusion of 
the higher order terms in the chiral expansion. Because the coefficients in the Lagrangian 
for QChPT are not well constrained, it is difficult to make numerical predictions of the size 
of the nonanalytic terms. However, it is certain that lattice simulations should be modified 
to account for this leading nonanalytic behavior. 
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